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In the scenario with Z mediated flavor changing neutral current occurring at the tree level due to 
the addition of a vector-like isosinglet down-type quark a" to the SM particle spectrum, we perform 
a \ 2 fit using the flavor physics data and obtain the best fit value along with errors of the tree level 
Zhs coupling, U s b- The fit indicates that the new physics coupling is constrained to be small: we 
obtain \U s t\ < 3.40 x 10 -4 at 3<r. Still this does allow for the possibility of new physics signals in 
some of the observables such as semileptonic CP asymmetry in B s decays. 
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I. INTRODUCTION 

The Standard Model (SM) of the electroweak interactions successfully explains most of the experimental data to 
date. However in recent years, there have been quite a few measurements of quantities in B decays which differ from 
the predictions of the SM. For example, inB^ ttK, the SM has some difficulty in accounting for all the experimental 
measurements The measured indirect (mixing- induced) CP asymmetry in some b — > s penguin decays is found 
not to be identical to that in B — > J/ipK s 0-0] , counter to the expectations of the SM. The measurement of indirect 
CP asymmetry in B s — > J / ijxf) by the CDF and D0 collaborations shows a deviation from the SM prediction 0-0] 1 ■ 
The observation of the anomalous dimuon charge asymmetry by the D0 collaboration [9|-|TT| also points towards some 
new physics in B s mixing that affects the lifetime difference and mixing phase involved therein [lzl . Il3j . A further hint 
of new physics has been seen in the exclusive semileptonic decay B — > K* u + ix~ : the forward-backward asymmetry 
(Afb) has been found to deviate somewhat from the predictions of the SM [14l - fl7j 2 . Though the disagreements are 
only at the level of ~ 2-3<r, and hence not statistically significant, they are intriguing since they all appear in b — > s 
transitions. Therefore the study of new physics effects in various b — > s observables is crucially important. 

A minimal extension of SM can be obtained by adding a vector-like isosinglet up-type or down-type quark to the 
SM particle spectrum [l9l - [34j |. Such exotic fermions can appear in E e grand unified theories as well in models with 
large extra dimensions. Here we consider the extension of SM by adding a vector like down- type quark d! . The 
ordinary Q em = — 1/3 quarks mix with the d' . Because the d' L has a different I3L from di,, sl and &l, Z- mediated 
FCNC's (ZFCNC) appear at tree level in the left-handed sector. In particular, a Zbs coupling can be generated: 

Zfcnc = - \ U sb &fP L bZ» + h.c. (1) 

Z COS uy/ 

This coupling leads to a new physics contribution to b — > s transition (such as B s -B s mixing, b — » sfi + /i~ & b — > svv 
decays, etc) at the tree level. This tree level coupling U s b can be constrained by various measurements in the b — > ,s 
sector. 

In this paper we consider observables such as B s -B s mixing, branching ratios of B — ¥ X s fi + fi~ , B s — > (J. + fJ,~ and 
B — > X s vv to constrain the new physics coupling U s b- Instead of obtaining the usual scatter plot which shows the 
allowed ranges of the U s b parameter space, we perform a x 2 fit which provides us the best fit value of U s b along 
with the errors. We then study the effect of tree level Zbs coupling on the indirect CP asymmetry in B s — > ipcf>, 
anomalous dimuon charge asymmetry a s sl , forward-backward (FB) asymmetry in B — > X s fj, + fi~ and the branching 
ratio of B s — > t + t~ . We show that the various measurements in the b — > s sector put strong constraint on the allowed 
values of U s b- However it is still possible to have new physics signals in some b — > s observables. 

The paper is organized as follows. In Sec. UH we discuss the methodology for the fit. In Sec. IIIIl we present the 
results of the fit. In Sec. IIVI we obtain predictions for various b —> s observables. Finally in Sec. [Vj we present our 
conclusions. 

II. METHOD 

As U s b denotes the Zbs coupling generated in the ZFCNC model, the parameters of the model are therefore the 
magnitude and the phase of this coupling, |i7 s b| and 4> s b = aigU s b- 

In order to obtain constraints on the new physics coupling U s b, we perform a x 2 fit using the CERN minimization 
code MINUIT [35}. The fit includes observables that have relatively small hadronic uncertainties: (i) the branching 
ratio of B — > X s /j, + fi~ in the low- and high-q 2 regions, (ii) the branching ratio of B s — > /i + /J~, (iii) the ratio of 
the branching ratio of B s — > and the mass difference in B s system, (iv) the branching ratio of B — > X s vv. 

We include both experimental errors and theoretical uncertainties in the fit. In the following subsections, we discuss 
various observables used as a constraint. 



1 The recent LHCb update does not confirm this result Q]. Their measurement is consistent with the SM prediction 

2 The recent LHCb update does not confirm this result [lSll . Their measurement of the Apg distribution is consistent with the SM 
prediction, except in the high-g 2 region. 
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B -> X s fi+ fx~ 



The effective Hamiltonian for the quark-level transition b — > s /i + /j in the SM can be written as 

10 



v i= 1 



(2) 



where the form of the operators Oi and the expressions for calculating the coefficients Ci are given in Rcf . [36] . The 
operator Oi, i = 1,6 can contribute indirectly to b — > s fi + fi~ and their effects are included in the effective Wilson 
coefficients Cg and C-j [H, [37[ . 

The Zbs coupling generated in the ZFCNC model changes the values of the Wilson coefficients Cg^o- The Wilson 



coefficients Cg°f in the ZFCNC model can be written as 



Ctot 
10 



C f n U sb (A „.J2 



C 



10 



oi V* s V tb 

7T U s b 



(4 sin^ 6 W - 1) 



(3) 
(4) 



Here V t * s V tb ~ -0.0403 e"* 1 . We use the SM Wilson coefficients as given in Ref. [37] 
The calculation of branching ratio gives 



BTZ(B -> X sl i+ pT) = 



a 2 BlZ(B -> V c e^) I^VJ 



tb\ 



A-K 2 f{rn c )n{m c ) |V cb | 5 



D{z)dz , 



where 



D(z) = (1 - z) 2 (1 + 2z) (|C* ot | 2 + |C^| 2 ) +4 1 + - \Cf\ 2 + 12Re(Cf C* ot *) 



(5) 



(6) 



Here z = q 2 /m 2 = (p^+ +P^-) 2 /w 2 and m q = m q /mb for all quarks </. The expressions for the phase-space factor 
/(m c ) and the 1-loop QCD correction factor ft(m c ) are given in (38l |. 

The theoretical prediction for the branching ratio of B — > X s /i + fi~~ in the intermediate q 2 region (7 GeV 2 < q 2 < 
12 GeV 2 ) is rather uncertain due to the nearby charmed resonances. The predictions are relatively cleaner in the 
low-g 2 (1 GeV 2 < q 2 < 6 GeV 2 ) and the high-g 2 (14.4 GeV 2 < q 2 < m 2 ) regions. We therefore consider both low-g 2 
high-g 2 region in the fit. 



We define \ 2 as 



where 



Anw = 



-Dhigh = 



(i- 



D{z)dz 
D(z)dz 



-^X s fi+ /A - :low 

BK(B ->■ X sf i+ fj,- 
BK(B -> X sf i+ ii 



f D\ ow 
V l 



5.69947\ 2 



82522 
£>hi B h- l-56735 \2 
0.635465 / 



J low 



4tt 2 f(m c ) K (m c ) \V cb \ 2 
a 2 BTl(B X D ev) \V*V tb \< 



5.69947 ± 1.82522 



(7) 
(8) 

(9) 



) high 2 



4tt 2 f(m c ) «(m c ) |V c6 | 2 



a 2 BK(B -> X c ev) \V t * s V tb \ 2 



1.56735 ± 0.635465 . (10) 



Here we have added an overall corrections of 30% to the theoretical prediction of B1Z(B — > X s /j, + fi )hi g h, which 
includes the non-perturbative corrections. 



B 3 



The purely leptonic decay B s — > /i + /i is chirally suppressed within the SM. The SM prediction for the branching 
ratio is (3.35 ± 0.32) x 10 -9 [48|. Recently LHCb collaboration reported a very strong upper bound on the branching 
ratio of B s -> p+ pT , which is 3.8 x 10~ 9 at 90% C.L. 



r/B = 0.5765 ±0.0065 [39] 
f bs = 0.229 ± 0.006 GeV [41, 42J 

Bi, s = 1.291 ±0.043 [41, 42] 
A M s = (17.69 ± 0.08) [45J 
'^ffi' 1 = 0.967 ± 0.009 [41] 
I V t * Vis | = -(0.0403 ± 0.0009) 
B1l{B -> XJLv) = (10.61 ± 0.17) x 10^ 



B1Z(B 3 -> fj, + n~) = (0.0 ± 2.30) x 10~ u [40] 
BTZ(B ^ X sf i+ fi-) low = (1.60 ±0.50) x 10 -6 [43,44] 
23ft(B X,jl* + /Owgh = (0-44 ± 0.12) x 10 -6 [43, 4^ 
Bft(B -» X s vu) = (0.0 ± 40) x 10" 5 [46] 
mt(mt) = 163.5 GeV 
m c /m b = 0.29 ± 0.02 
t Bs = (1.520 ± 0.020) ps' 1 [8J 



TABLE I: Inputs that we use in order to constrain \U a b\-^ab parameter space, when not explicitly stated, we take 
from Particle Data Group [47^ . 



The branching ratio of B s — > n + ji in the ZFCNC model is given by 



We define x 2 as 



with 



BK{B S -> = 



GWM B mlfl s r B 
16tt 3 



v; s v tb \\ l- 



4m 2 , 



M I ^rtot 1 2 



M% 



|C^ t | 2 -0. 0\2 
13.5408 



\C 



tot 1 2 
10 I 



16w 3 BK{B S ->■ (i+fj,-) 



G 2 F c?M Bs mlflT Bs \V* s V tb \\l - 



0.0 ± 13.5408 



Mi 



C. Ratio of BTL{B S -¥ m + /0 and the mass difference in the B s system 

The mass difference AM S is given by 

AM S = 2|Aff 2 M | . 

The SM contribution to Mf 2 is 

r s,SM _ /ta*t/ \ 2 i\*2 



M 



12 _ ^W* V ^ 2M w MB ^Bf B B BB E(xt) , 



where x< = mf/M^y and ?7b is the QCD correction. The loop function E(xt) is given by 



E{x t ) 



-Axt + 11x1 - x t 3xjlnx t 



4(1 -x t Y 2(1 -x t ) 3 

The mass difference AM S in the ZFCNC model is given by [28[ 



A W - = 0l^^| 2 M^M s ^B/ b 2 s i3 6s |i?(^)||A s 



A s is given by 



where 



A=i + „|jLl-if ^ b 



.<7(a; t ) , 2V2tt 2 
a = 4-^ — r , o ■ 



E(xt) ' G F M 2 / S(a; t ) 
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The loop function C(xt) is given by [28 



C(x t ) 



4 



3xt In xt 



(20) 



l-z t (l-x t ) 2 _ 

The term in Eq. jT^I proportional to a is obtained from a diagram with both SM and new physics Z vertices; that 
proportional to b corresponds to the diagram with two new physics Z vertices. 
Dividing Eq. (HD by Eq. ([17]), we get 



We define x 2 as 



^M^ B B bs \E{x t ) 



Xbr- 



ly^t Ot I 2 

l°10 I 

|A 3 | 



0.0 N 



13.6328 



4m 2 |Ctot|2 
Ml \A S \ 



(21) 



(22) 



with 



AM,. 1 



4m J 
Ml" 



0.0 ± 13.6328 



(23) 



D. B -> X.i/S 



The effective Hamiltonian for the decay B X s vv is given by 



H eff = y t ',Mofe)Wv-A(Hv-A + h.c. , (24) 

V2 27r sin £V 



with 



Xt 



2 + x t 3xt — 6 



8 ixt — i (x t - ly j 

The presence of tree level Zbs coupling changes the value of the structure function X (x t ). The structure function 
within the 2TCNC model can be written as 

X^x t )=X (x t )+(^^)u sb . (26) 
V aV ts Vtb J 

The branching ratio of B — > X s vv is given by [4!| 

*^ w=B ^ M snls)' (27) 

where C 2 is given by 



C = o , • 4. \V t * s V tb X>(x t )\ 2 . (28) 



27r 2 sin 



We define x 2 as 



with 



/ \Vt*Vt b X>(xt)\ 2 -0.0 \ 2 
Xb^Xsuu y 0.069157 J ' 1 j 

IV* V V'/ ^l 2 2 7 r 2 sin 4 ^|^| 2 /(m c )«(m e ) 
B1Z(B — ¥ X c ev) ryar 

Here we have used the present upper bound B7Z(B —> X s vv) < 64 x 10 -5 at 90% C.L. (4(| which can be written as 
(0.0 ±40) x 10~ 5 . 

Therefore the total \ 2 can be written as 

Xtotal = Xg^x s fj,+ /i-:low + X§^X S fj,+ ju- :high + ^B t ->/i+/j- + XBR-mix + Xb^X s uv ' 
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Parameter Value 

\U^\ (0.90 ± 0.83) x 10" 4 

<j) sb (0.00 ± 181.34)° 

x'/d.o.f. 1.72/3 



TABLE II: The results of the fit to the parameters of ZFCNC model. 



Observables 


Predictions 


SM 


ZFCNC 


(j>f (rad) 





(0.00 ± 0.03) 


|A S | 


1 


1.01 ±0.01 


<i x 10 s 


(1.92 ±0.67) 


(1.98 ± 13.88) 


Br(B a ->• r+r+) x 10 Y 


5.74 ±0.27 


3.34 ± 1.92 


(g^r GeV^ 


3.33 ±0.25 


3.38 ± 0.26 



TABLE III: ZFCNC predictions for potential observables. 



III. RESULTS OF THE FIT 



The results of these fits are presented in Table HU It may be observed that the x 2 per degree of freedom is small, 
indicating that the fit is good. We observe that the present flavor data put strong constraint on Zbs coupling. At 3er, 
we obtain \U sb \ < 3.40 x 1CT 4 . 



IV. PREDICTIONS 



A. Semileptonic asymmetry a s sl 



The expression for the semileptonic asymmetry a s sl is given by 



IP? 



M- 



12 



\M[f M \ |A fl | ' 



(32) 



where the CP violating phase 4> s is defined by the following equation, 



Arg 



M 



12 



(33) 



The parameter A s takes into account the new physics effects in mixing and is defined as 



M{ 2 =M^ SM (1 



M 



M 



.NP 

^) = Mt 2 SM A s = M? 2 SM \A s \e<> 

12 



(34) 



Thus 6° can be written as 



(35) 



where 



ASM 



(3.84 ± 1.05) x 10" 3 [51J. Also, one has [52, M 



|r? 2 



\M*f M \ 



(5.0 ± 1.1) x 10 



-3 



(36) 



The predictions for 4>f, |A S | and a s sl in ZFCNC model are given in Table HID We see that it is possible to have large 
deviations in <f> s (and hence a s sl ) from its SM predictions. 
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B. Zero of Forward-Backward asymmetry 

The FB asymmetry of muons in B — > X s /i + fj,~ is obtained by integrating the double differential branching ratio 
( , d d gTC a ) with respect to the angular variable cos 9 [HH 



f 1 dcosO , d2 , BU , - f , d cos 61- 

JO a cos y J — 1 c 



a / \ Jij az a cos y j— i dcosO 

fndcose-^Bn+f.dcose-j^ ' 1 j 

JO «2 dcosv J— 1 (12 a cost? 

where # is the angle between the momentum of the B-meson and that of /i + in the dimuon center-of-mass frame. 
Within the ZFCNC model, FB asymmetry in B — > A s /j + is given by 

^fb(z) = , (38) 

where D(z) is given in Eq. |6]and E(z) by 

£(z) = Re(C^ ot C 1 t g t *)^ + 2Re{CfC[f*) . (39) 

Zero of Afb(z) is determined by 

E{z) = Re(C* ot C*^ t *)z + 2Re(Cf r C 1 ° t *) = . (40) 



The prediction for (g 2 ) ncl in ZFCNC model is given in Table Hill One can see that large deviations from SM 
prediction is not possible. 



C. BK{B S -+ t+ O 
The branching ratio of B s — > t + t~ in the ZFCNC model is given by 



3a 2 T Bs m 2 T / 4m 2 . |Cj; ot ' 2 

8wM^B^\E{^t)\ Y 1 ~ M% ~~\A S 



BK(B S -> r+r") = _„, 2 L ,^ 1 - T7^^AM S (41) 



The prediction for B1Z(B S — > t + t ) in ZFCNC model is given in Table Hill We see that it is possible to have large 
suppression in B7Z(B S — > r + r~) as compared to its SM prediction. 



V. CONCLUSION 



In this paper, we consider a minimal extension of the SM by adding a vector-like isosinglet down-type quark d! 
to the SM particle spectrum. As a consequence, Z-mediated FCNC's appear at tree level in the left-handed sector. 
In particular, we are interested in Zbs coupling which leads to a new physics contribution to b — > s transition such 
as B s -B s mixing, b s/i + /i~, b — > svv decays, etc at the tree level. Using inputs from several observables in 
flavor physics, we perform a \ 2 n t to constrain the tree level Zbs coupling, [7 sb . The fit takes into account both the 
theoretical as well as the experimental uncertainties. 

We conclude the following: 

• x 2 P er degree of freedom is small, indicating that the fit is good. This is expected as the SM itself is in good 
agreement with the data. 

• The present data put strong constraint on the Zbs coupling. At 3cr, \U s b\ < 3.40 x 10~ 4 . 

• Despite the strong constraint on the Zbs coupling, it is possible to have new physics signals in some b — > s 
observables such as semileptonic CP asymmetry in B s decays. 
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